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Abstract 



Frequency distributions of local muon densities in high-energy extensive air-showers 
(EAS) are presented as signature of the primary cosmic ray energy spectrum in the 
knee region. Together with the gross shower variables like shower core position, angle 
of incidence, and the shower sizes, the KASCADE experiment is able to measure 
local muon densities for two different muon energy thresholds. The spectra have 
been reconstructed for various core distances, as well as for particular subsamples, 
classified on the basis of the shower size ratio N^/N^. The measured density spectra 
of the total sample exhibit clear kinks reflecting the knee of the primary energy 
spectrum. While relatively sharp changes of the slopes are observed in the spectrum 
of EAS with small values of the shower size ratio, no such feature is detected at 
EAS of large N^/Nq ratio in the energy range of 1-10 PeV. Comparing the spectra 
for various thresholds and core distances with detailed Monte Carlo simulations the 
validity of EAS simulations is discussed. 

Key words: cosmic rays; air shower; muon component, energy spectrum, mass 
composition 
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1 Introduction 



Measurements of the energy spectrum and the elemental composition of the pri- 
mary cosmic radiation constrain theoretical models of the sources, acceleration 
mechanisms and transport of the radiation through the interstellar space. While 
for lower energies direct measurements by satellites or balloon-borne detectors 
yield spectroscopic results (see ref.[l]), for primary energies above some 10^^ eV 
only indirect measurements via extensive air shower (EAS) observations can be 
performed. It is well known that the energy spectrum of the primary cosmic ra- 
diation shows a kink (mostly referred to as "knee") at energies around 3 PeV [2]. 
Though the first evidence of the existence of this knee has been presented more 
than 40 years ago [3], the knowledge of the detailed structure of the spectrum 
in the PeV region is still scarce, and the origin of the knee not yet understood. 
Most of the earth-bound air shower experiments use large detector arrays to 
measure charged particles and reconstruct shower sizes of the individual events 
by adjusting a particular lateral distribution function to the measured densities. 
The resulting shower size spectra reflect the primary energy spectrum, but a 
quantitative conversion to energy has to invoke a model of the shower develop- 
ment and on an assumption of a mass composition. Hence the determination of 
the energy spectrum is affected by different systematic uncertainties, especially 
by the dependence on the model of high-energy interactions. This also leads to a 
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mutual dependence of the results for the energy spectrum and mass composition. 
When comparing recent results of earth-bound air shower experiments, like 
CASA-MIA [4], TIBET [5] or Akeno [6], significant differences in the absolute 
magnitudes of the total flux, the knee position and slope of the energy spectrum 
are noticed. A recent non-parametric analysis of KASCADE data [7] reports 
equally large differences on the energy spectra, depending on the high-energy 
hadronic interaction model, illustrating the considerable influence of the inter- 
action models underlying the Monte Carlo simulations. To identify the ultimate 
sources of the disagreements, it would be useful to analyze different experiments 
on basis of a coherent methodology as well as to compare the resulting features 
for various sets of different EAS parameters in the individual experiments. 
In the present paper we endeavor to analyze the frequency distribution of local 
muon densities at fixed distances from the shower core. The local muon density 
spectra reflect the gross features of the primary energy spectrum, as the muon 
content for a certain distance to the shower center observed at sea-level is mainly 
determined by the primary energy. While the reconstruction of electron or muon 
size spectra necessarily implies a choice of the form of the lateral distribution 
function, spectra of the muon density are free from this bias. Thus with "inde- 
pendent" measurements of such spectra for different fixed core distances allow a 
check on the lateral distribution obtained from simulations. In addition, the lay- 
out of the KASCADE experiment [8], with a central detector system consisting 
of densely packed muon counters with different shielding, enables the study of 
density spectra for two different muon energy thresholds. Hence the consistency 
of the simulations with respect to the muon energy spectrum can be performed. 



2 Experimental setup and data handling 

KASCADE (KArlsruhe Shower Core and Array DEtector) is a multi-detector 
setup [8] at Forschungszentrum Karlsruhe (110m a.s.L), Germany, for EAS mea- 
surements in the primary energy range around the knee. The main detector 
components of KASCADE used for the present analysis are an "array" of 252 
stations, located on a squared grid with 13 m spacing and a "central detector" 
comprising additional detector systems. 

The array is organized in 16 subarrays (4x4 stations each) and provides the data 
necessary for the reconstruction of the basic EAS characteristics like electron and 
muon size (total number of electrons and muons in the EAS), core location, and 
arrival direction of individual air showers. The special arrangement of shielded 
and unshielded detectors on top of each other allows an independent estimation 
of the total electron and muon number for each individual shower. The densities 
are estimated and corrected iteratively for punch-through effects (/x-counters) 
and muon contamination (e/7-counters). Lateral correction functions from sim- 
ulations of EAS and detectors are used. The reconstructed particle densities 
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are fitted by Nishimura-Kamata-Greisen (NKG) functions in tlic experimental 
accessible distance ranges (10-200 m for the electron component, 40-200 m for 
the muon component). The densities obtained are then integrated from zero to 
infinity for the total numbers and A^^. In addition we quote the so-called trun- 
cated numbers for which the NKG functions are integrated in a limited range 
only {Nf = Pe ■ 2nRdR and N^^ = J^o^^ Pn ■ ^nRdR). These truncated 

numbers provide reduced systematic uncertainties since extrapolations into the 
radial range outside our measurement areas are avoided. Uncertainties are esti- 
mated by Monte Carlo calculations and range below 20% for the total numbers 
and even better for the truncated numbers. The location of the shower core once 
inside a fiducial area is determined to better than 3 m. The arrival direction of 
the shower is reconstructed from the arrival times of EAS particles (cr ~ 0.5°). 
These procedures are described in detail elsewhere [11]. 

The KASCADE central detector is placed at the geometrical center of the de- 
tector array. It consists of four different detector systems (Fig. 1), covering a 
total area of 16x20 m^. The local muon density of EAS is measured with the 
multiwire proportional chambers (MWPC) and the trigger plane. 
Below the hadron calorimeter [9], with a total thickness of 5 cm lead, 154 cm 
steel, and 77 cm concrete corresponding to a threshold for vertical muons of 
2.4 GeV, a setup of 32 large multiwire proportional chambers is installed [10]. 
The chambers of three different sizes (5m^, 8.2 m^, and 8.8 m^) are arranged in 
two layers with 38 cm vertical distance. In total each layer has a sensitive area 
of 129 m^. A single chamber consists of two layers of cathode strips at angles 
of ±34° with respect to a layer of anode wires in between. The chambers are 
operated with an argon-methane gas mixture. The electronics allows a digital 
readout of all wires and strips, i.e. in total of c. 32,000 channels. Hits of through- 
going particles are reconstructed as the intersection of the anode wires and the 
cathode strips. Thresholds and delays arc adjustable for each channel separately. 
A continuous monitoring of the reconstruction efficiency during measurements 
is performed. The chambers have a spatial resolution of about 5 mm. 
The reconstruction of high-energy muons in the MWPC starts from the recon- 
structed hits in each plane and the shower direction. The reconstructed direction 
is required to agree with the shower direction within ±15° in zenith and ±45° in 
azimuth (the azimuth cut is not used for showers with zenith angles of < 10°). 
These cuts appear reasonable as for core distances below 100 m high-energy 
muon tracks are nearly parallel to the shower axis. It is known from simulations 
that in the considered range of primary energy (PeV) and core distances the 
muon density (for > 2.4 GeV) very rarely exceeds 1 per m^. Therefore recon- 
struction ambiguities are negligible. High-energy 5-electrons which are produced 
to a small amount in the absorber are eliminated by calculating the height of the 
intersection of two nearby tracks. If they cross inside the central detector, the 
track with the larger deviation from the shower axis is rejected while the other 
is accepted as a muon only. About 0.1% of tracks are rejected by this cut. The 
spatial resolution of single tracks is about 1.0 cm, the angular resolution is ~ 1°. 
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Fig. 1. Schematic view of the KASCADE central detector. It consists of four detector 
systems: an 8-layer hadron calorimeter [9], scintillation counters at the trigger plane 
and top cluster, and a setup of multiwire proportional chambers (MWPC) [10] below 
the calorimeter. 

The number of tracked muons N* is also corrected for the reconstruction effi- 
ciency which is estimated for each single data acquisition run 12 h) separately. 
The efficiency was found to be very stable with a mean value of (e) = 93% [10]. 
The local muon density p* for each EAS is defined by N* divided by the total 
sensitive area A* of the MWPC setup. Due to the layout of the chambers, A* 
depends on the angle of incidence of the shower and is calculated for each event 
individually {{A*) = 107 m^ for the selected EAS). Only that area where muons 
parallel to the shower axis would penetrate the whole absorber and both cham- 
ber planes, is taken into account for the calculation of the muon density. 
The second detector system is a layer of 456 plastic scintillation detectors in 
the third gap of the calorimeter, called trigger plane [9]. Each detector consists 
of two square plates of plastic scintillators (47.5x47.5x3 cm^) separated by a 
wavelength-shifter, which is read out by a single photomultiplier. Fast electronics 
records low-energy (muons) and high-energy deposits (cascading hadrons) and 
provides a trigger for the calorimeter and other detector systems. In the present 
analysis the trigger plane with an active area of 208 m^ is used to estimate the 
local density of muons with a threshold of 490 MeV for vertical incidence. 
The muon density p|f at the trigger plane is reconstructed in the following way: 
To remove signals from cascading hadrons in the absorber an upper limit of the 
energy deposit of 30 MeV in each of the 456 scintillation counters is imposed. 
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Detectors with larger energy deposits and their immediate neighbours are not 
considered for further reconstruction. For the remaining detectors, the energy de- 
posit and the sensitive area, both corrected for the shower direction, are summed 
up. The number of reconstructed muons is then calculated by the sum of 
the energy deposits divided by the mean energy deposit of a single muon in the 
shower. According to Monte Carlo calculations this mean value depends slightly 
on the core distance Rc ((-Edcp) = 7.6211 — 0.00495 ■ i?c in MeV) and is corrected 
for. Rc is the distance in meter of the core position to the center of the trigger 
plane (or MWPC) projected to a plane perpendicular to the shower axis. The 
density p|f is obtained as ratio of and the sensitive area of the trigger plane 
for each individual event {{A^^) = 202 m^). 

The MWPC setup is triggered by trigger plane and top cluster, but not by the 
detector array. But the array, trigger plane and the top cluster are triggering all 
other components. The trigger plane fires if more than 7 detectors have signals 
> 1/3 mip or if at least one detector has E^^p > 300 MeV (for single hadron 
detection). The top cluster triggers if more than 8 (out of 50) detectors show a 
signal. An array trigger is activated if half of the stations of at least one subar- 
ray show an energy deposit (> 1/3 mip). One of the central detector triggers in 
conjunction with the array trigger have to be active to initiate the event recon- 
struction. 

After some general cuts (core position less than 91 m from array center, < 40°, 
lg{Nj^) > 4.745 — 0.212 Ig(Ai'e)), more than two million events have been used 
for the present analysis, recorded in circa 282 days of measuring time. Mea- 
sured fluxes have been corrected by 9% for the dead-time of the data acquisition 
system. 



3 Local muon density spectra 

3.1 All-particle spectra 

It is reasonable to assume that at a fixed distance from the shower axis the lo- 
cal muon densities map the energy of the primary particles [11] and that muon 
density spectra carry information about the primary energy spectra. 
The reconstruction of muon density spectra have been performed for two energy 
thresholds and for nine core distance ranges (Figures 2 and 3). These ranges 
are chosen in such a way that the sampling area are of equal size (1473.4 m^) 
and large enough to get reasonable statistical accuracy but retain small system- 
atic uncertainties due to the extension of the core distance bins. To suppress 
punch-through effects of the hadronic or electromagnetic component, EAS with 
Rc < 30 m are excluded. EAS with Rc > 72 m are excluded, too, because 
they can have their core outside the KASCADE array {Ro < 91m) if they are 
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Fig. 2. Integral spectra of the local muon density p* as measured by the MWPC system 
for different core distances Rc (from top to bottom with increasing Rc). The upper 
limits of the radial bins are 37.0, 42.9, 48.0, 52.7, 57.0, 60.9, 64.7, 68.2, and 71.6 m, 
respectively. The lines represent the results of the fit procedure (see text). The dashed 
line displays for one case the result of a fit with a single power law. 

very inclined. The p^-spectra are affected by trigger efficiencies at low densi- 
ties. The limitation at high densities are given by reconstruction uncertainties. 
For the MWPC system (higher energy threshold) these uncertainties begin at a 
fixed muon density of p* = 0.6 m~^ due to ambiguities in the track reconstruc- 
tion and punch-through effects of cascading hadrons. The hadronic energy and 
hadron particle density in EAS are increasing similarly to the muon density at 
all core distances [11]. In case of the trigger plane the cut on the deposited en- 
ergy (30 MeV) in each scintillation detector has a systematic influence on local 
densities above p^^ = 0.8 m~^. Here the density will be reduced since the intrinsic 
density fluctuations in the EAS together with fluctuations in the energy deposit 
lead to detector signals exceeding 30 MeV without punch-through effects. 
Figure 2 and Figure 3 show the flux spectra for the two muon thresholds in inte- 
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Fig. 3. Same as Figure 2 but for the local muon density /o)? measured by the trigger 
plane. Here the dashed lines display, for one case, the fit function of a two power law 
fit (below and above the knee position) without an intersecting knee region. 

gral form. The flux values are multiplied by (p*)^'^ and (pj^'')^", respectively. All 
spectra show a slight, but significant kink with decreasing density for increas- 
ing core distance. For the fit procedure the flux lg{j^) is assumed to follow a 
power law below and above the knee region. The following form of the differential 
spectra is assumed: 



hi 



h + PMp,) for /(7(p«) < lg{p,) < lg{p^^)) 
a[b-lg{p,)f + c for lg{p^^^) < lg{p,) < lg{pf) 
b2 + f32lg{p,) for /^?(p(f)) < lg{p^) < IgipW) 



The fit procedure estimates the indices (3i of these power laws, the position of 
the knee (if existing), and the boundaries of the different regions. The values of 
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the boundaries are estimated by the method of finite and dividing differences 
[12]. Especially the lower thresholds of the spectra, where they begin to deviate 
from a power law dependence are defined by this method. With increasing core 
distances these trigger thresholds move to lower muon densities because of the 
decreasing lateral distribution function. The position of the knee is calculated as 
the weighted center of gravity of the bins inside the knee region. The fit functions 
are included in Figures 2 and 3, as well as the particular position of the bending. 
The "width" of the knee region for all spectra amounts to /S.lg{p^/m~'^) ^ 0.15 . 




Fig. 4. Knee positions of the muon density spectra vs. core distance for both energy 
thresholds. The error bars indicate the uncertainties of the fit procedure. 

It is a remarkable observation that within the statistical uncertainty the knee 
positions for all spectra occur at the same EAS fiux, for both thresholds and 
all core distances (at a differential flux of d'N/dlg{pfj,/m~^^) = (1.2 ± 0.2) ■ 10~® 
m~^s~^sr~^ ). Figure 4 displays the knee positions for all reconstructed spectra. 
If we assume that the knee is a feature of the primary energy spectrum, the data 
points mark the lateral distribution of muons for two different energy thresholds 
for a flxed primary cosmic ray energy. This shows that, for the EAS registered, 
the fraction of muons between 490 MeV and 2.4 GeV is small. Figure 5 shows the 
power law indices of the density spectra. A higher muon energy threshold results 
in steeper spectra. This indicates a comparatively larger increase of the muon 
density per primary energy interval with increasing muon energy threshold. The 
spectra for different core distances are almost parallel leading to nearly constant 
indices for a given muon energy threshold. This conflrms previous experimental 
results [11] of only slight changes of the shape of the muon lateral distributions 
with increasing primary energy (which is different for the electromagnetic com- 
ponent of EAS). For both energy thresholds there is a clear difference in the 
indices below and above the knee. 
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Fig. 5. Power law indices of the differential muon density spectra for both energy 
thresholds. The error bars indicate the uncertainties of the fit procedure. 

Several tests were performed to check the robustness of the shape of the recon- 
structed spectra against experimental uncertainties or systematic features of the 
analysis. Differently chosen functional forms for the fit procedure to the differ- 
ential spectra result in similar indices of the power laws. The significance of the 
knee remains stable (see also dashed lines in Figure 3). The assumption of a 
single power law for the spectra (e.g. dashed line Figure 2) leads to reduced x^- 
values of around 1.5 — 2., whereas the used procedures fit the data with x^-values 
close to one. In addition, it has been found that effects of the binning of the den- 
sity and of the core ranges, or of the azimuthal distribution of the showers are 
negligible and within the statistical uncertainty of the spectra. Variations of the 
chosen zenith angle range shift the total spectrum in slightly, but the form of 
the spectra remains stable. This is reasonable, because differently inclined show- 
ers (of the same primary energy) generate slightly different numbers of muons, 
while the variation of the density with energy does not change. 

3.2 Spectra of EAS subsamples 

The ratio of the muon to electron content of EAS is traditionally considered as 
a mass-sensitive observable [13,14], since heavy ion induced EAS tend to have 
a large ratio due to the faster development of the electromagnetic component 
in the atmosphere. The shower sizes N^, and A'"*'' have been reconstructed for 
each individual event from the data of the array stations and the EAS have 
been classified according to the ratio and then divided in "electron-rich" and 
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"electron-poor" . The shower sizes are converted to the sizes of vertical showers 
to eliminate the influence of the different zenith angles: 

ln{N'J = ln{Ne) - f- ■ (sec ^ - 1) 
ln{N-) =ln{Np-^^^. {sec e-1) 

where A^e, Nj^, and 9 are the reconstructed quantities of the EAS, and X — 
1022 g/cm^ is the observation level. The quantities Ag and A^ denote the ab- 
sorption lengths of the electron and muon components in the atmosphere. The 
values were obtained from Monte Carlo simulations and parameterised as Ag = 
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Fig. 6. Examples for measured spectra of different muon content. The "all" -particle 
spectra have already been shown in Figures 2 and 3 and are here compared with the 
spectra of "electron-poor" and "electron-rich" EAS for the same core distance range. 
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104.3 + 13.5 ■ lg{Ne) g/cYv? and = 5 ■ Ae . Especially the electron number 
depends significantly on the zenith angle due to the rapidly increasing atmo- 
spheric absorption. The separation of the total sample of EAS in "electron- 
rich" and "electron-poor" showers is performed by a cut in the ratio l^ratio = 
lg{N'^) /lg{N2) = 0.75 . This value is optimized by Monte Carlo calculations (see 
section 4.1). The classification is done independently using the local muon den- 
sities at the central detector. For both subsamples the spectra are deduced in 
the same way as the "all-particle" spectra. 

As example Figure 6 shows the reconstructed local muon density spectra for 
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Fig. 7. Power law indices of the muon density spectra for both energy thresholds 
and all core distance ranges for the samples of electron-rich {f3i below and above the 
knee) and electron-poor (Ph) EAS. The error bars indicate the uncertainties of the fit 
procedure. 



(-Rc) = 45.5m for all, for the electron-rich (predominantly light ion induced), 
and for the electron-poor (predominantly heavy ion induced) showers. The gen- 
eral features of the spectra are similar for all core distance ranges; the compo- 
nent of electron-rich EAS dominates the flux below the knee while it strongly 
decreases after the kink. No knee is seen in the component of electron-poor EAS, 
and the spectra can be described by a single power law. The resulting slopes of 
the spectra, especially the differences of the slope-values for the two thresholds 
and subsamples, are very similar for the various core distances as shown in Fig- 
ure 7. Whereas the assumed fit functions describe the all-particle spectra well, 
the spectra for the electron-rich EAS are not well described by power laws above 
the knee. Also for the electron-poor sample slight deviations from a pure power 
law dependence at the high energy end are observed. This holds for all radial 
ranges. An energy dependent separation efficiency of the primary masses as well 
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as astrophysical sources (composition, acceleration, propagation) can cause these 
deviations from simple power law dependencies. 



4 Comparisons with simulations 

4.1 Air-shower simulations 

For the interpretation of the measured muon density spectra in terms of the 
primary energy spectrum a-priori knowledge inferred from Monte Carlo simula- 
tions of the air-shower development is necessary. The present analysis is based 
on CORSIKA simulations (version 5.62) [15] and a full and detailed simulation 
of the detector response. The simulations have been performed using the inter- 
action model QGSJET [16] for the high-energy interactions and GHEISHA [17] 
for interactions below E'lab — 80 GcV and subsequent decays. The electromag- 
netic part of the showers is treated by EGS4 [18]. Observation level, earth's 
magnetic field, and the particle thresholds are chosen in accordance with the 
experimental situation of KASCADE. The U.S. standard atmosphere [15] was 
adopted. The simulations cover the energy range of 5 • 10^^ - 3.06 • 10^^ eV di- 
vided into 7 overlapping energy bins with a spectral index of —2.7 . For each bin 
200 showers are simulated except for the two highest energy ranges where only 
100 and 50 showers were generated, respectively. The calculations are performed 
for three zenith angular ranges (0° — 15°, 15° — 20°, 20° — 40°) and for three 
primary masses: protons, oxygen and iron nuclei. The response of KASCADE is 
simulated by a detailed detector simulation program based on the GEANT [19] 
package. Each generated shower is passed ten times through the detector sim- 
ulation. The shower cores are randomly distributed over the KASCADE array 
within a circular area of 95 m radius around the center. Hence, a total statistics 
of 103,500 EAS is used. The output of the simulations is analyzed by the same 
procedures as applied to the measured data, reducing systematic uncertainties. 
Figure 8 displays examples of function of Eq for different muon thresh- 

olds, core distances, and primary masses. The selection cuts have been applied 
as to the measured data. The error bars indicate the width of the distributions. 
They decrease with increasing energy and mass. A power law dependence is fit- 
ted in a restricted energy range to reduce the influence of showers with primary 
energies outside the simulation range. 

Figure 9 shows the resulting power law indices S for all core distances in case of 
primary protons and iron nuclei. The slopes are nearly independent of the ra- 
dial distance, confirming the weak dependence of the shape of the muon lateral 
distribution with primary energy [11]. But the slopes depend on muon threshold 
and primary mass. Systematic uncertainties are at the 10% level and indicated 
in Figure 9. These systematics are estimated by varying the energy and angular 
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Fig. 8. Examples of local muon densities vs. primary energy of simulated EAS for 
different muon thresholds, core distances, and primary masses. The error bars indicate 
the standard deviations of the densities. The lines show power law fits taking into 
account the statistical uncertainty of the mean values which are smaller than the 
marker sizes. 
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Fig. 9. Power law indices 5 of the cc Eq relation for different Rc and muon thresholds 
in case for primary protons and iron nuclei. The error bars indicate the systematic 
uncertainty (see text). 

spectrum of the simulated events, as well as modifying the degree of fluctua- 
tion of the observables resulting from the simulations. The uncertainty of the 
high-energy interaction model itself is tested with a set of simulations based 
on the VENUS high-energy interaction model [20]. The slope of the oc Eq- 
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Fig. 10. Distributions of the parameter l^atio used for mass separation shown for 
measured (left hand) and simulated (right hand) EAS. In case of simulations the 
energy dependence of the ratio is shown. The error bars represent the r.m.s. of the 
distributions. The dashed lines indicate the dividing line between "electron-rich" and 
"electron-poor" . All selection cuts are applied. 

dependence differs systematically by ~ +0.1 for both energy thresholds, but this 
systematics is not included in the error bars in Figure 9. 

The detailed simulations allow also to verify that the cut on the shower size 
ratio yratio = h{^'ij,)/^9{^e) is energy independent. The distribution of Yy-aHo 
for the measured events, and the quality of the mass separation as provided by 
the simulation calculations are shown in Figure 10. To divide the total sample 
in electron-rich and electron-poor EAS a cut of Fratio = 0.75 is chosen. By this 
most of the proton induced showers belong to the "electron-rich" class, whereas 
primary iron and medium nuclei are associated to the "electron-poor" class. Ob- 
viously the classification is nearly energy independent (Figure 10 right). Also, 
the qualitative behavior of the electron-poor and electron-rich distributions was 
found to be insensitive to small changes of l^atio ■ 



4-2 Features of the energy spectra 



When relating the density spectra to the primary energy spectrum of cosmic 
rays a power law spectrum ^ oc Eq'^ is assumed. The energy spectrum can be 

written as ^ • where ^ has to be deduced from the EAS simulations and 



dpf, dEo ' dEo 



(X (pfj) ^ is taken from the experimental results. Thus the spectral index 7 



can be expressed by ^ — S ■ {P — 1) + 1 with 6 from the simulations (p^ oc E(f , see 
Figure 9). If the correct elemental composition is adopted, all measured muon 
density spectra (of the total sample or of a certain subsample) should result 
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consistently in the true primary energy spectrum, irrespective which core dis- 
tance and muon energy threshold are considered. Hence by use of the results 
of various core distance ranges and different muon energy thresholds systematic 
effects induced by the Monte Carlo simulations could be checked. E. g., a pos- 
sible dependence of the slopes and of the knee position on core distance would 




core distance R^. [m] 



Fig. 11. Variation of the reconstructed power law index below the knee of the pri- 
mary all-particle energy spectrum with the core distance for both energy thresholds 
assuming a pure proton composition. 

indicate an insufficient description of the muon lateral distribution by the sim- 
ulations while a comparison of the spectra observed with different muon energy 
thresholds could reveal inconsistencies of the simulated muon energy spectrum. 
As an example, Figure 11 shows the variation of the resulting exponents 71 of 
the primary all-particle spectrum below the knee derived for both muon energy 
thresholds under the assumption that the primaries are protons. The smaller 
fluctuations in with increasing atomic number result in a decrease of the 
spectral index of the resulting primary energy spectrum for the assumption of 
a heavier composition. Figure 12 compares the knee positions determined under 
the assumptions of proton primaries with the results found for the case assum- 
ing iron primaries. The knee positions resulting from iron nuclei as primaries are 
systematically shifted to smaller energies as compared to a pure proton compo- 
sition. This is due to the fact that the local muon density is increasing with the 
primary mass (see Figure 8 left). 

The density spectra for the different core distances are independent of each other 
and the resulting slopes and knee positions of the primary energy spectrum agree 
within their statistical uncertainties. This supports the confidence in the lateral 
distribution predicted by the Monte Carlo simulations, and allows to present 
results averaged over all core distance bins (Table 1). Nevertheless there remain 
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Fig. 12. Variation of the estimated knee position in the primary all-particle spectrum 
with the core distance for both energy thresholds, with the assumption of a pure 
proton or pure iron cosmic ray beam. 

obvious systematic differences in the results for the two muon energy thresh- 
olds, observed for all core distances. The systematic differences might arise from 
possibly incorrect assumptions on mass composition due to the sensitivity of 
the muon spectrum to primary mass. Such an effect, however, should be consid- 
erably reduced when analysing the electron-rich and electron-poor subsamples 
which should be enriched in light and heavy primaries, respectively. For these 
samples no variation with core distance is again observed, and Table 1 presents 
average values. But the systematic differences for the two thresholds remain. 
In order to check the influence of heavier contributions the electron-rich EAS 
sample has been additionally analyzed assuming larger fractions of helium and 
even oxygen primaries. A flattening of the spectra and a shift of the knee po- 
sition to lower energies occurs by up to 5% (with the extreme assumption of 
100% oxygen). When assuming a pure oxygen composition for the electron-poor 
sample the primary spectrum steepens by A7 ?a 0.15. Such effects do not ex- 
plain the systematic discrepancy displayed by the results from the two different 
muon energy thresholds. Therefore we conclude that an incorrect description of 
the muon energy spectrum by the Monte Carlo simulations is the origin of the 
discrepancy. 

The effect does not only occur for the QGSJet model used for the present analy- 
sis. A smaller sample of reference showers generated with the VENUS model has 
been used to study the observed difference. A general shift to a steeper primary 
energy spectrum (A7 ^ 0.2) and a lower knee position is found. That may be 
associated to differences in the modelling of the high-energy interaction [7] . How- 
ever, the inconsistency with respect to the two different muon energy thresholds 
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Table 1 

The spectral indices, fluxes and positions of the knee in the primary energy spectrum 
for the different EAS samples analyzed assuming single element primaries {dJ]inee/dE 
is given in [m~^s~^sr~^GeV~^]). The parameters do not depend on the core distance. 
Therefore mean values are given. 



all-particle spectrum 




for primary proton 


for primary oxygen 


for primary iron 


it 
ll 

dJLJdE 


2.68 ±0.02 
2.89 ±0.02 
6.61 ±0.05 
(10.4 ±2.0) -10-1^ 


2.66 ±0.02 
2.85 ±0.01 
6.53 ±0.05 
(12.5 ±2.5) -lO-i^ 


2.51 ±0.02 
2.71 ±0.01 
6.47 ±0.05 
(13.9 ±2.9) -10-14 


tp 
7i 

tp 

72 

lg«ce/GeV) 
dJljdE 


2.89 ±0.01 
3.11 ±0.02 
6.48 ±0.05 
(3.0 ±0.8) -10-14 


2.96 ±0.01 
3.18 ±0.01 
6.43 ±0.04 
(2.3 ±0.3) -10-14 


2.75 ±0.01 
2.96 ±0.01 
6.37 ±0.04 
(3.7 ±1.1) -10-14 





electron-rich sample 


electron-poor sample 




for primary proton 


for primary iron 


ll 

72* 

lg(^knee/GeV) 


2.83 ±0.03 
3.41 ±0.03 
6.70 ±0.05 


2.40 ±0.01 


tp 

7i 
tp 

72 

lg(4L/GeV) 


2.97 ±0.02 
3.53 ±0.02 
6.53 ±0.04 


2.69 ±0.01 



persists. The considered muon energies are comparatively low, and are treated 
in the CORSIKA simulation code mainly by the low-energy interaction model 
GHEISHA. Thus the inconsistencies are most probably due to the low-energy 
model. There are in fact indications for deficiencies of the code from another 
study [21]. 



5 Summary and conclusions 



Frequency spectra of local muon densities of EAS in the PeV region were mea- 
sured and analyzed for various core distances and for two muon energy thresh- 
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olds. For both thresholds the all-particle spectra show the knee structure, i.e. two 
power laws with increasing steepness in the knee region. Compared to shower 
size spectra based on the charged particle or electron number [22,23], the muon 
density spectra show a relatively smooth knee with a small, but clear change of 
the power law exponent. 

With help of the muon to electron number ratio, estimated on an event-by-event 
basis, the registered EAS are divided into electron-rich and electron-poor sub- 
samples. The subsample of the electron-rich EAS shows the same knee features 
as the total sample but with a more pronounced knee. The electron-poor sample 
shows no change of slope within the density range investigated. An identical 
feature has also been observed in combined energy and composition analyses 
of size spectra measured by KASCADE [24] albeit then based on Monte Carlo 
simulations. Simulations indicate that electron-rich showers originate from light 
primary nuclei. This feature holds irrespectively of details of the interaction 
models. Hence we conclude that the knee reflects a feature of the light particle 
spectrum and that the spectrum of heavy particles does not change slope in the 
range of our measurements. Such a behaviour is expected if the knee is caused 
either by interstellar magnetic fields or a change of the interaction in the atmo- 
sphere since in these cases the knee of nuclei should be displaced by a factor of 
Z or A, respectively, to that of protons. 

Detailed EAS and detector simulations were used to interpret the measured 
muon density spectra in terms of the primary energy spectrum. Independent 
of the elemental composition assumed all measured spectra should result in the 
same primary energy spectrum, irrespectively of core distance or muon threshold. 
However, only when assuming the true composition the derived energy spec- 
trum will be the correct one. This agreement for different spectra actually is 
observed for varying core distances. Hence we conclude that the muon lateral 
distribution is sufficiently well represented by the simulations. In contrast, the 
results for the two energy thresholds lead to different exponents and positions 
of the kinks. These differences are larger than the systematic uncertainties due 
to the unknown composition, especially in the case of the subsamples. Thus, 
the measurements presented here reveal that the Monte Carlo simulations are 
not capable to describe the muon energy distribution correctly. Such indications 
arise also from studies of the muon lateral distributions for different muon en- 
ergy thresholds with KASCADE [25]. In view of these systematic discrepancies, 
it is difficult to draw definite conclusions but some general features of the pri- 
mary energy spectrum can be stated: The all-particle energy spectrum exhibits 
a knee at -Ekncc ~ (3 — 5) ■ 10^^ eV with a change of the spectral index of order 
A7 0.2 — 0.3 . This knee is only seen in the light ion subsample, at the same 
position but with a distinctly larger steepening of A7 fa 0.5 . The heavy ion com- 
ponent of the cosmic ray flux displays no steepening in the energy range of 1-10 
PeV and a smaller slope than the light component below the knee. Within the 
uncertainties the findings about the all-particle spectrum are compatible with 
recent results from KASCADE [7] and other experiments [4-6]. 
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From the experimental point of view the study of muon density spectra es- 
tabhshes a new approach to investigate the energy spectrum of cosmic rays. 
Although the statistical accuracy is limited the measured spectra reflect the fea- 
tures of the primary energy spectrum in an astonishingly direct manner. But con- 
siderable inconsistencies arise when attempting to convert the measured muon 
density spectra into the primary energy spectrum based on simulations. It is 
only due to multiparameter measurements of experiments such as KASCADE 
that the deficiencies of the simulations get revealed. Thus improvements of the 
hadronic interaction models incorporated in the simulations appear to be the 
most important prerequisite for a consistent interpretation of the data in terms 
of elemental composition and energy spectrum of primary cosmic rays in the 
knee region. 
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